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Abstract 

We analyze the structure of networks minimizing the global resistance to flow (or dissipated energy) with 
respect to two different constraints: fixed total channel volume and fixed total channel surface area. First, we 
determine the shape of channels in such optimal networks and show that they must be straight with uniform 
cross-sectional areas. Then, we establish a relation between the cross-sectional areas of adjoining channels at 
each junction. Indeed, this relation is a generalization of Murray's law, originally established in the context 
of local optimization. Moreover, we establish a relation between angles and cross-sectional areas of adjoining 
channels at each junction, which can be represented as a vectorial force balance equation, where the force 
weight depends on the channel cross-sectional area. A scaling law between the minimal resistance value and 
the total volume or surface area value is also derived from the analysis. Furthermore, we show that no more 
than three or four channels meet in one junction of optimal bi-dimensional networks, depending on the flow 
profile (e.g.: Poiseuille-like or plug- like) and the considered constraint (fixed volume or surface area). In 
particular, we show that sources are directly connected to wells, without intermediate junctions, for minimal 
resistance networks preserving the total channel volume in case of plug flow regime. Finally, all these results 
arc illustrated with a simple example, and compared with the structure of natural networks. 



Networked structures arise in a wide array of different contexts such as water, gas and power supply 
of a city, vascular systems of plants and animals, or river basins PQ[2][II|- Thus, optimization of transport 
in networks has evident industrial and economical importance, but may also shed light on the structure 
of natural networked structures. Indeed, the analysis of these structures from optimization and selection 
principles has been recently the subject of intense scientific activity @] |5] [S] |E| and controversy |§||1(J | [TT ] . 
Besides, theoretical models - based on local optimization (i.e. optimization of the geometry of a single 
junction) - have been attempted to explain in detail the regular patterns of vascular networks |12 j |13 | [T3 ) |15 j . 
However, it is generally known that as the global optimum is achieved the local optimum of a single junction 
is often discarded. In the present paper, we characterize the structure of networks satisfying to the global 
optimization of transport. For the class of networks mentioned here, euclidean metric must be taken account, 
and the optimization must be achieved with respect to some geometrical constraint. 

Precisely, the problem we consider can be expressed as it follows: consider s sources at the same potential 
(electrical potential, pressure, concentration, temperature,...) Vs and w wells at the same potential Vw, their 
respective positions being fixed. What is the architecture of the network linking all the sources to all the 
wells and minimizing the effective resistance (or dissipated energy), for a fixed total channel volume or fixed 
total channel surface area ? Or equivalently, which architecture minimizes the total channel volume or 
surface area for a same value of the global resistance ? 

In the following, we shall refer often to the electrical circuit terminology, although this study obviously 
concerns any flow-in- network situation. Let us denote each pipe by a pair of indices corresponding 
to the labels of its two ends. We suppose a priori that pipes can be curved, but we assume that their 
aspect ratios are sufficiently high so a length 1^ and a local cross-sectional area s«(/) (where I denotes the 
curvilinear coordinate along a channel) can be unequivocally defined for each pipe The resistance dr^ 



I 



of an infinitesimal piece of pipe of length dl is then defined as: 

dtin = —dl, (1) 

where p is the "resistivity", supposed to be the same for all the pipes. For m = 1, the flow in each channel 
is plug-like, while for m — 2 the flow is Poiseuille-like. Assuming there is no leakage through the pipe lateral 
surface, the resistance of the whole pipe is: 

hi 

Tn = f ^dl. (2) 

Since we shall inspect the minimal resistance configuration with respect with two different constraints (a 

fixed total channel volume V tot and a fixed total surface-area channel Stot), we introduce for simplicity the 

hi 

"constraint function": C n = J sfjdl , so that: C\ — Vtot, and Ci/ 2 oc Stot- 

(j,i)o 

1 Cohn's theorem 

To characterize the architecture of minimal resistance networks, we shall invoke Cohn's theorem, originally 
developed in the context of electrical circuit analysis |17) : consider a one-port network composed entirely 
of two-terminal elements with resistances r%j . The variation of the effective network resistance R with the 



variation of the resistance ry is given by 
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(3) 



No particular assumption is made on the expression of the resistances r*y for the derivation of this result 
(indeed, the theorem is still valid for complex impedances). Conservation of flow and energy only are required. 
Thus, Cohn's theorem can be applied to a broader class of flow-in-network situations. 
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2 Optimal shape of channels 

We first notice that in order for the effective network resistance to be at its minimum value with respect to the 
constraint C„, each channel must be straight with a uniform cross-sectional area (i.e.: Sij(l) — s,j). Indeed, 
we see from eq. [SJthat the effective network resistance R is a monotone function of the individual resistances 
Tij. Thus, any small change in pipe diameter or pipe length from the minimal resistance configuration - 
compatible with the constraint - must lead to an increase of the resistances r^ . As a consequence, the 
length of each pipe must be as small as possible and its diameter as large as possible, i.e. each channel must 
be straight with a uniform cross-sectional area. Besides, it can be noticed that a circular cross-sectional 
area have the specific property of minimizing both the pipe surface area for a fixed volume (or equivalently 
maximizing the pipe volume for a fixed surface area) and the dissipative energy in the channel for a fixed 
incoming flow-rate in case of Poiseuillc-flow regime. 

3 Relations between diameters: generalized Murray's law 

We now establish relations between diameters and angles in an optimal network, for a given topology (meaning 
that no junction or channel can be added or removed from the network, but the channel lengths and cross- 
section areas are free to vary). We thus have to minimize the function R = R + AC„ (where A is a Lagrange 
multiplier) with respect to the independent variables {sij} and {r^ = (a;,, j/i)}, respectively the channel cross- 
sectional areas and node positions. Using Cohn's theorem |3| the condition of extremum with respect to the 
cross-sectional areas (dR/dstj = 0) gives: 




(4) 
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Furthermore, conservation of flow-rate at each junction i (/Jiy = 0) implies: 



(m+n)/2 



= 0. 



This relation, illustrated on Fig. ^ 
Murray's law ^2] to any flow profile 
the particular case m = 2, n = 1). IN 
optimization of the network structu 
optimization (flow and channel crossj 
was found for a given flow, and not 
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Figure f : Relation between cross-sectional areas of adjoining channels in a minimal resistance network. This 
relation is a generalization of Murray's law. 



4 Geometry of nodes 

Condition of extremum with respect to the node positions (dR/dri = 0) together with relation ^straightfor- 
wardly leads to the following vectorial equality at each node i: 



where is the outward-pointing unit vector along the channel (i,j) (see Fig. |2J). This equality, relating 
angles between adjoining channels to their cross-sectional areas, is similar to a force balance equation, where 
the weight of the force acting along the channel is directly proportional to s™- . As for Murray's law, local 
optimization principles have already been proposed in order to describe the geometry of nodes in natural 
networks, namely: minimization of channel volume (V), channel surface area (S), dissipated power (P), and 
drag force (D) on the walls (HEUd]. ^ these approaches consist in varying the position of a given junction, 
while the positions of the other junctions, the network topology, the channel cross-sectional areas, and the 
flow-rates through every channel arc remained fixed. However, in the context of a global optimization, a 
change in a node position should alter the flow-rate distribution, and it is therefore to be expected that 
global minimization of the dissipated energy leads to a different optimal geometry of nodes than in the local 
optimization context (P). Indeed, the optimal geometry of nodes described by Eq. HJis similar to the one 
obtained for (S) (when n = 1/2) or (V) (when n = 1), but different from (P) |g||T5 | |Tl | . 

5 Scaling-law between minimal resistance and constraint value 

A relation between the minimal resistance value and the constraint value can be established, using Eq. 21 
and conservation of energy: 




(6) 



3 




(7) 



5 



Figure 2: Relation between angles and cross-sectional areas of adjoining channels in a minimal resistance 
network. This relation is similar to a force balance equation describing the equilibrium of strings tied together 
and under respective tensions, or weights, sfj. 

On the other hand, a classical result of optimization theory relates the Lagrange multiplier to the change 
of the minimal resistance with respect to the constraint value: A = — ^^7 (note that Eq. implies A > 0). 
Therefore, it is found that the resistance of an optimal network scales as C„ m ^ n , i.e.: 



where I is a parameter with dimension of length, depending solely on the network topology, the positions of 
sources and wells, and the values of m and n. 

We have shown that a minimal resistance configuration, for a given topology, if it does exist, must satisfy 
to the equations 01 El and 03 Wether the extrema characterized by this set of equations are local minima 
or local maxima is not clear (although this uncertainty might be dispelled by some convexity argument). 
Nevertheless, because individual resistances have finite values, there must exist at least one configuration 
with global minimal resistance (but we do not know if this configuration is unique) |18| . 




(8) 



6 Upper bound on the node connectivity 

Finally, we establish an upper bound on the number of channels joining in one node, in a bi-dimensional 
minimal resistance network. To do so, we look at a given junction of N channels and determine when this 
junction is preferentially replaced with two junctions respectively of 3 and N — 1 channels. Suppose we create 
a new channel of infinitesimal length dl^, as depicted in Fig. |3J Then, the length variation of the two other 
channels joining in the new 3-fold junction are: dl\ = —dl^cosOi and dli = —dl^ cos92, with: 6\ + 82 = 7, 
where 7 is the angle between these two adjacent channels. The variation of the associated resistances are 
respectively: dr\ = — pdl$ cos#i/s™, d- r i — ~ pdh cos^/s™, and dr^ — pdl^/s™, where si, S2 and S3 are the 
respective channel cross-sectional areas. Moreover, this transformation must preserve the value of C„, so the 
new channel cross-sectional area S3 must satisfy: 



Suppose now that the iV-fold junction was in a minimal resistance configuration. Then, conditions 01 and 
must be fulfilled, and we can replace i\ and i\ in Eq. EUby their expressions (Eq. Moreover, conservation 
of flow rate relates ^3 to i\ and 12'- «3 = — ii — 12- Using Eq. we see that the resistance variation dR is 
negative when: s ^ m+n )/ 2 > s ( m +")/ 2 -|- s ' m+ ™^ 2 . The sign in the right-hand side of this inequality is positive 
when the two adjacent channels are crossed by flows in same direction, and negative when they are crossed by 
flows in opposite directions. The former inequality can be rewritten as: cos#i + r n cos #2 > ( 1 i r~ ^~ j + , 
with r = S2/S1. Before establishing an upper-bound on the node connectivity, we must notice the following 
"rules" on the geometry of junctions in optimal networks: 



s™ cos 9 1 + s 2 l cos 9 2 • 



(9) 



Using once again Cohn's theorem, we obtain the variation of the effective resistance: 




(10) 
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Figure 3: Elementary transformation of a TV-fold junction to a (N — l)-fold junction plus a 3-fold junction. 
A new channel, with infinitesimal length dl% is thus created. 

1. There is at least one angle lower than 2-k/N between two adjacent channels in a iV-fold junction (from 
geometrical consideration) . 

2. There is at least one pair of adjacent channels crossed by flows in opposed directions (from flow con- 
servation) . 

3. The angle between two adjacent channels is always lower than ir (from Eq. |SJ). 

Let us choose 9\ and O2 such that sin#i = r"sin#2, what corresponds to the maximum value of the 
left-hand side of the former equality. Since 7 < ir (rule 3.), we easily check that both 6\ and 62 are positive 
and lower than 7r/2, and simple algebra leads to: 

cos6»i + r" cos 6» 2 = \J\ + r 2n + 2r n cos 7. (11) 

Thus, the resistance variation dR is negative if and only if: 

(l ± r :R P L ) m+n -l-r 2n 
cos 7 > f ± (r) = ± '— , (f 2) 

where the functions / + (r) and /_ (r) correspond to the respective situations of two adjacent channels crossed 
by flows in same and opposite directions. The analysis of /+ (r) and /_ (r) shows that, for any value of r, 
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these functions are bounded in the following way: /+ (r) < 2 m + n — 1 for any value of m, n; /_ (r) < if 



m > n, and /_ (r) < — 1 if m = n. So if 7 is lower than 7+ = arccos ^2 m +» — 1 j for the first situation, or 

7_ = 90° (if m > n) or 180° (if m = n) for the second situation, we are ensured that the resistance variation 
is negative. Let us inspect the different situations: 

• If m = 2, and n = 1/2: 7+ ~ 97.4°, 7_ = 90°. We know there is at least one angle lower than 
360°/TV between two adjacent channels in a TV-fold junction (rule 1.). By choosing this angle as 7 in 
the previous analysis, we conclude that a TV-fold junction is preferably replaced with a (TV — l)-fold 
junction plus a 3-fold junction, as long as TV > 4. The new structure is not in a minimal resistance 
configuration, determined by Eqs0]and[Sl so the "relaxation" of the new structure such a configuration 
implies a further decrease of the effective resistance. Eventually, we can repeat the same reasoning on 
the (TV — l)-fold junction, if TV — 1 > 4. We come to the conclusion that exactly three channels meet 
at each junction in such an optimal network. 

• If m = 2, and n = 1 or m = 1, and n = 1/2: 7+ ~ 74.9°, 7_ = 90°. Following the same argumentation, 
we conclude that a TV- fold junction is preferably replaced with a (N — l)-fold junction plus a 3-fold 
junction as long as N > 5. Thus, no more than four channels meet in one junction in such an optimized 
network. Furthermore, it can be noticed that only two kinds of 4-fold junctions can exist in such a 
network: either three adjacent channels are crossed by flows of same sign (and the last flow has an 
opposed sign), or two adjacent channels are crossed by flows with same sign and the two other adjacent 
channels are crossed by flows with same opposite sign. A 4-fold junction with channels crossed by flows 
with alternate signs is preferably replaced by two 3-fold junctions, since there is always two adjacent 
channels crossed by flows with opposite signs and with an angle lower than 90° (rule 1.). 

• If m — 1 and n = 1: 7+ = 0°, 7_ = 180°. But we know that there is always two adjacent channels 
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crossed by flows with opposite signs in a AT- fold junction (rule 2.), with an angle between them lower 
than 180° (rule 3.). So the TV-fold junction is preferably replaced with a (N — l)-fold junction plus a 
3- fold junction for any N > 4. Now, if we let the new structure of the network "relax" to a minimal 
resistance configuration, it must simultaneously satisfy Eqs. [5]and[B]at every junction, and particularly 
at the 3-fold junction. But this set of equations applied in a 3-fold junction has only trivial solutions 
when m < n: either one cross-section is null, or the three channels are colinear. We conclude that 
sources are directly connected to the wells, with no intermediate junction, in a minimal resistance 
network preserving total channel volume and in case of plug-flow regime. 

As a concluding remark for this section, we point out that the same reasoning may be used on the total 
channel length variation instead of resistance variation (Steiner tree problem). In that case, we obtain that 
links meet at threefold junctions (with equal angles of 120°) in a length-minimizing network. 

7 A simple example 

We compare our results with a simple example: two sources and two wells placed at the corner of a rectangle, 
as depicted on Fig. 01 Four configurations are analyzed. In configuration (1), sources are directly connected 
to wells, without any intermediate junction. In Configuration (2), sources are connected to wells via a 
4- fold intermediate junction. In configurations (3) and (4), sources and wells are connected through two 
3-fold junctions (and the position of intermediate junctions are chosen such that equality [fj] is satisfied). The 
corresponding dimensionless resistance (~^) " °f each configuration is reported on tabled From these 
expressions, we note the following observations, in agreement with our results: firstly, we notice that R scales 
as (1/C n ) m . Secondly, when m — n (= 1), configuration (1) is the smallest resistance configuration, for 
any value of the aspect ratio b/a. Thirdly, resistance of configuration (3) is always lower than resistance of 
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configuration (2) and higher than resistance of configuration (1) (Ri < R3 < R2), for any value of to, n, 
and b/a. Fourthly, resistance of configuration (4) is lower than resistance of configuration (2) as soon as: 
\J 1 + (b/a) 2 < 2 2 (" l ~™)/(" l +") ; for any value of to, n, and b/a. One can easily check that this criterion on the 
aspect ratio b/a (for given values of to and ri) corresponds to the condition for Eq. [5] to be simultaneously 
satisfied with Eq. 0at each 3- fold junctions of configuration (4). In particular, resistance of configuration 



(4) cannot be lower than res: 
Fifthly, when m > n, resista] 
sufficiently low value of b/a. 




ffl (2) 




Figure 4: Four different network configurations linking two sources to two wells, placed at the corners of a 
rectangle of length a and b. 
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Configuration 


pa \ a / 


(1) 




(2) 


2^(l + (!) 2 ) " 


(3) 


/ — 






(4) 


(4 + t 2 ^)"(ur 2 (i) f ), 




with li = -* ^ — 



Table 1: Dimensionless resistances J| ( ^ ) " corresponding to the four configurations depicted on Fig. 0] 
For configurations (3) and (4), l\ = l\ja and l<i = h/a are the dimensionless lengths of the two kind of 
channels. 

8 Comparison with natural networks 

All the results derived in previous sections (relations QJ [S] H3 |H1 as well as the upper-bound on the node 
connectivity) are consequences of global optimization. However, these results have been established by 
studying any local perturbation of the structure. Such a local adaptive process may take place during 
ontogeny of natural networks. Therefore, it may be of interest to compare the structure of some natural 
networks with the results presented in this work. Indeed, it has been already shown in various publications 
|19 j |2U | that Murray's law is well satisfied in some appropriate portions of human and animal vascular 
systems. In that case, the flow profile is nearly Poiseuille-like (m = 2) and the relevant constraint is 
a fixed total channel volume (n = 1) ^HHOI- Validity of Murray's law for vascular system of plants is 
more controversial [20| 123 E3 E3I > mostly because of the underlying theoretical assumptions in the original 
derivation of Murray's law, and of the particular structure of veins in vascular system of plants. Nevertheless, 
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experimental data suggest that a relation N sign (iy) = is still verified, with v between 2.49 and 3 

3 

|15j |23 | [21 | |251. Let us look more precisely at the bi-dimensional leaf venation network, like the one reported 
on Fig. |SJ Leaf veins are actually vascular bundles |26j . supporting two parallel flows : a pressure-driven flow 
of water and minerals from petiole to stomata through xylem tissues, and a diffusive flow of nutrients and 
photosynthesis products in the opposite direction through phloem tissues. So petiole (or major vein) and 
stomata play alternatively roles of sources and wells for the leaf. An outer layer of cells, called the bundle 
sheath, surrounds the vascular tissues. Although this layer is not fully impermeable, the leaky radial flow is 
small when compared with the axial flow, except for the minor veins |21 j |27 ) . For these veins, leakage is very 
important and the pressure field and nutrient concentration nearby are almost uniform. Kull and Herbig 28 
investigated on leaf topology of several species. They observed that leaf venations preferably show trivalent 
nodes with six neighbors, and noticed that this geometry is typical of self-generating structures like bubble 
floats. In a recent study, S. Bohn et al. |29| analyzed geometry of junctions in the leaf venation of various 
species. They observed that angles between veins are very well defined and that a vectorial balance equation 
comparable to eq. can be established, where the weight of each vector is directly proportional to the vein 
radius (i.e. n = 1/2). Comparison of Bohn et al. observations with our optimization principles suggest then 
that structure of leaf venation corresponds to a minimization of the resistance for a fixed total channel surface 
area (or minimization of surface area for a fixed value of resistance) . This result is coherent with the idea 
of a predominant building cost of the bundle sheath cells over those of the vascular tissues [2] [2H] ■ Taking 
n = 1/2 and comparing Eq. with experimental studies of Murray's law leads to a value of to between 
1.99 and 2.5, meaning that flow in veins is nearly Poiseuille-like. Note it is assumed in the theory that all 
channels have same resistivity p. However, density of xylem and phloem tissues in a leaf vein might be a 
function of the vein diameter as well. This variation of the resistivity is then include in the coefficient to, 
what could explain the slight difference observed between the experimental value of to and the theoretical 
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value for a Poiseuillc-fiow regime. All these results and the presence of trivalent nodes suggest then that 
structure of leaf venation correspond to minimal resistance configuration preserving the total surface area 
and nearly Poiscuillc flow profile. The measure of the scaling-law between the hydraulic resistance and the 
total channel volume or surface area might be an additional way to test this conjecture. 




Figure 5: Portion of leaf venation. In most species, the structure is netted-like, and veins meet in 3-fold 
junctions. 
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